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Abstract 

In this paper I deal with the WDO 137-349 binary system consisting 
of a white dwarf (WD) and a brown dwarf (BD) in a close circular 
orbit of about 116 niin. I, first, constrain the admissible range of values 
for the inclination i by noting that, from looking for deviations from 
the third Kepler law, the quadrupole mass moment Q would assume 
unlikely large values, incompatible with zero at more than 1-sigma level 
for j ^ 35 deg and ? > 43 deg. Then, by conservatively assuming that 
the most likely values for i are those that prevent such an anomalous 
behavior of Q, i.e. those for which the third Kepler law is an adequate 
modeling of the orbital period, I obtain z = 39 ± 2 deg. Such a result 
is incompatible with the value i = 35 deg quoted in literature by more 
than 2 sigma. Conversely, it is shown that the white dwarf's mass 
range obtained from spectroscopic measurements is compatible with 
my experimental range, but not for i = 35 deg. As a consequence, my 
estimate of i yields an orbital separation of a = (0.59 ± O.O5)R0 and 
an equilibrium temperature of BD of Teq = (2087 ± 154)K which differ 
by 10% and 4%, respectively, from the corresponding values for z = 35 
deg. 
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1 Introduction 

In general, binary systems composed by a white dwarf (WD) orbited by 
a brown dwarf (BD) as a companion are rare: the WD0137-349 (BPS CS 
29504-0036) system (Maxted et al., 2006), consisting of a WD and a BD 
orbiting in a 116 min close circular path, belongs to such a class. BD must 



have survived a previous phase in which it was engulfed by the red giant 
progenitor of WD (Pohtano, 2004) experiencing orbital drag which notably 
shrank its orbit, originally much wider. 

In (lorio, 2007) a dynamical determination of its quadrupole mass mo- 
ment Q from deviations of the third Kepler law was claimed obtaining an 
unlikely large value for it: Q = — 1.5 x 10^'' kg m^, with \Q/MR^\ 10^, 
where M and R are the WD's mass and equatorial radius. In this paper I 
clarify the meaning of such a result and show how it can be used to con- 
strain the inclination angle i. My analysis furnishes a physical justification 
of the use of the third Kepler law in modeling the orbital period of this 
particular system, without assuming it uncritically a priori. Moreover, it 
also corrects an error concerning i in (Maxted et al., 2006) in which i = 35 
deg is quoted, yielding an orbital separation of O.65R0. Such values for i 
and a are reported in other works being used for investigations on various 
aspects of the WD0137-349 system like, e.g., the determination of Q itself 
(lorio, 2007) and the heating of BD (Burleigh et al., 2006). 

2 Constraining the inclination with the quadrupole 
mass moment 

The phenomenologically measured period amounts to (Maxted et al., 2006) 



in principle, it accounts for all the dynamical effects affecting our binary 
system to the measurement accuracy. Let us, now, contrast it to the purely 
Keplerian period 



where a is the relative semimajor axis, and M and m are the masses of WD 
and BD, respectively, in order to see if it is compatible with Pb) within the 
errors, or if some significant discrepancy occurs. Such a comparison does, 
indeed, make sense because the "ingredients" entering eq. (2) have all been 
determined in a way which is independent of the third Kepler law itself, 
apart from the inclination i which will, thus, be treated as a free parameter. 
Indeed, the mass of WD (Maxted et al., 2006), 



Pb = 0.083 ± 0.0002 d : 



(1) 




M = (0.39 ± 0.035)M, 



(3) 
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was non-dynamically inferred from its effective temperature Tefj and surface 
gravity log (7 (Driebe et al., 1998; Benvenuto and Althaus, 1999)- measured, 
in turn, from an analysis of the hydrogen absorption lines of the optical spec- 
trum of WD (Koster D. et al., 2001)-using well tested and reliable models 
of white dwarfs. Then, from 

KmPh _ _ Ma sin i 
It: [M + m) 

where K^j^ are the projected semiamplitudes of the radial velocities, phe- 
nomenologically measured from the spectroscopic radial velocity curves, 
^M/m the projected barycentric semimajor axes and i is the inclina- 
tion of the orbital plane to the plane of the sky, it is possible to measure the 
ratio of the masses 7 = m/M as (Maxted et al., 2006) 

7 = = 0.134 ±0.006 : (5) 

I have used (Maxted et al., 2006) 

i^M = 27.9.5 ± 0.3 kms~\ 

(6) 

Km = 208.5 ±8.1 km s"^ 

From the value of M and eq. (5) the mass of the BD follows (Maxted et al., 
2006) 

m = (0.052 ±O.OO7)M0. (7) 

Now, from 

« = (1±7)^, (8) 
smz 

it is possible to express P^^^ in terms of known quantities, apart from the 
inclination i which will be treated as an independent variable on which I 
want to put some dynamical constrains. Note that in (Maxted et al., 2006) 
the value i = 35 deg is released yielding an orbital separation a = O.65R0 
which, in turn, was used to obtain an equilibrium temperature of BD of 
(Burleigh et al., 2006) 

jHBD) ^ ^(WD) ^ 2OOOK. (9) 
As we will see, such a value for i is probably wrong. 
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Let us investigate the ratio 

AP _ Pb-i^^°^ 

5{AP) ^ 6Pb + <5P(o) 

as a function of i. From Figure 1 it can be noted that outside the interval 
35 deg < i ^ 43 deg such a discrepancy becomes significative at more that 
1-sigma level changing its sign as well: inside such a range it is, instead, 
compatible with zero. 

May such a discrepancy have a physical meaning, so that the associated 
figures for i can be considered realistic values for it? Deviations from the 
third Kepler law could be induced by dynamical effects neglected in modeling 
the orbital period. However, post-Newtonian corrections of order 0{c~'^) 
(Soffel, 1989) are negligible because much smaller than the experimental 
error = 0.0002 d. Another possible candidate is the post-Keplerian, 
Newtonian effect of the quadrupole mass moment Q of WD (lorio, 2007). 
Let us explore this possibility in order to see if it is viable; such an analysis 
will also better clarify the strange result obtained in (lorio, 2007) for i = 35 
deg. We will see that it can yield us useful insights about the inclination of 
the considered binary system. 

By using eq. (23) of Appendix 

37r 3 

in Figure 2 I plot Q/6Q and Q/M versus where 5Q is given by eq. (25)- 
eq. (26) of Appendix. It can be noted that Q, assumed as responsible of the 
discrepancy between Pb and P'^'^\ becomes, in fact, a determined quantity 
at more than 1-sigma level outside about 35 deg ^ ^ ^ 43 deg: WD would 
be highly oblate {Q < 0) for « < 35 deg and highly prolate {Q > 0) for i > 43 
deg. The physical meaning of such an outcome is, however, suspect because 
Q/MR^, where R = (0.0186 ± O.OO12)R0 is the WD's radius (Maxted et 
al., 2006) would get as large as 150-200. Such large values of Q/MR^ are 
unlikely; WD's quadrupole moment (normalized to MR?) are expected to 
be smaller than 1 (Baym et al., 1971; Papoyan et al., 1971). In principle, it 
could be argued that the determined Q is, in fact, an "effective" quantity 
which accounts for other, unknown dynamical effects. A more conservative 
interpretation of our results is that the inclination of the orbital plane of the 
WD0137-349 system is confined in that interval in which the third Kepler 
law is an adequate description of the orbital dynamics and Q is compatible 
with zero. 
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Figure 1: Discrepancy AP between the measured orbital period and the 
computed Keplerian one P^^\ normahzed to its error 6{AP), as function 
of the inchnation i. It becomes significant, at more than 1-sigma level, for 
i < 35 deg and i ^ 43 deg. 
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Figure 2: The adimensional quantities Q/5Q (upper panel) and Q/ME? 
(lower panel) as functions of the inclination i. The upper curve in the lower 
panel is determined hy Q + 6Q, while the lower curve in the lower panel is 
determined by Q — 6Q, so that the admissible values for Q/MB? lie in the 
region delimited by such two curves. 
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I will, now, work within such a framework. Since I have information 
on M independent of the third Kepler law itself, I can use them to further 
constrain i. Indeed, by equating Pb to P^^^ and using eq. (4) and eq. (8) I 
have 

'-(l + 7mPb^'/'^ 



arcsin 



2ttGM 



39 ±2 deg. (12) 



Such a value is incompatible with z = 35 deg quoted in (Maxted et al., 2006) 
at more than 2-sigma level. The relative semimajor axis, computed with eq. 
(8) and eq. (12), becomes 

a = (O.59±O.O5)R0, (13) 

yielding an equilibrium temperature for BD of 

rj^^D) = (2087 ± 154)K. (14) 

It is interesting to see what are the dynamical constraints on M, obtained 
by modeling the orbital period with the third Kepler law, and compare them 
to the spectroscopic ones yielding eq. (3). Indeed, from eq. (2), eq. (4) and 
eq. (8) I get 

GM = {l + ^f(^){^)\ (15) 
\27r/ \smz / 

used to draw Figure 3 in which the region delimited by Af it 5M is displayed 
as a function of i. Note that the admitted range of mass values for i = 35 
deg is 0.4457 < M < 0.5757, which does not overlap with the spectroscopic 
one of eq. (3). 



3 Discussion and conclusions 

In this paper I got dynamical constraints on some orbital and physical pa- 
rameters of the WD0137-349A/B binary system consisting of a white dwarf 
and a brown dwarf as a companion. By using the correction to the third 
Kepler law due to the quadrupole mass moment Q, I was able to put con- 
straints to the inclination angle i of the system. By discarding those values 
of i which would yield significant deviations from the third Kepler law asso- 
ciated to unlikely large Q I obtained i = 39 it 2 deg. The spectroscopically 
inferred range of admissible values for the mass M of the white dwarf is 
compatible with our determination of i which also yields an orbital separa- 
tion of a = (0.59 lb 0.05)Rq and an equilibrium temperature of the brown 
dwarf of Teq = (2087 it 154)K. An incorrect value quoted in (Maxted et al., 
2006) of i = 35 deg yields, as a consequence, a = O.65R0 and Teq 2000K. 
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Figure 3: The allowed values for M (Solar masses), obtained by modeling 
the orbital period with the third Kepler law, lie in the region delimited by 
the upper curve (M + 6M) and the lower (M — 6M) curve. 
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Appendix: the contribution of the quadrupole mass 
moment to the orbital period 

One of the six Keplerian orbital elements in terms of which it is possible to 
parameterize the orbital motion in a binary system is the mean anomaly M 
defined as Ai = n{t — Tq), where n is the mean motion and Tq is the time of 
periastron passage. The mean motion n = 2tt/P\^ is inversely proportional 
to the time elapsed between two consecutive crossings of the periastron, 
i.e. the anomalistic period Pb- In Newtonian mechanics, for two point- 
like bodies, n reduces to the usual Keplerian expression rv- > = 
where a is the semi-major axis of the relative orbit and M = nii + m2 is 
the sum of the masses. In many binary systems, as WD0137-349, the period 
Pb is accurately determined in a phenomenological, model-independent way, 
so that it accounts for all the dynamical features of the system, not only 
those coming from the Newtonian point-like terms, within the measurement 
precision. 

Here I wish to calculate the contribution of the quadrupole mass moment 
Q to the orbital period in a more general and accurate way than done in 
(lorio, 2007) by retaining the orbital eccentricity. According to lorio (2007), 
the acceleration induced by the quadrupole mass moment Q can be cast into 
the form 

The quadrupole mass term Aq is small with respect to the usual Newto- 
nian monopole term Aq^ so that it can be treated perturbatively. In order to 
derive its impact on the orbital period Pb; let us consider the Gauss equa- 
tion for the variation of the mean anomaly in the case of an entirely radial 
disturbing acceleration A 

dM 2 , fr\ (1 -e^) ^ ^ 

-n A[-]+- -Acosf, (17) 



dt na \aJ nae 

where / is the true anomaly, reckoned from the periastron. After inserting 
Aq into the right-hand-side of eq. (17), it must be evaluated onto the 
unperturbed Keplerian ellipse 

^ (18) 



1 + e cos / 



By using (Roy, 2005) 

(^)'v^- (19) 



dM \r. 
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eq. (17) yields 

df n(l + ecos/)^ 
di " (1 - e2)3/2 



1 + 



3GQ(1 + ecos/)= 
71.^0^(1 — e^)'^ 



The orbital period can be obtained as 



(1 - e2)3/2 /-a- r 3GQ(l + ecos/)= 



n 



From eq. (21) it can be obtained 
with^ 



cos /(I + e cos /) 
2^ 



cos /(I + e cos /) 
2^ 



(20) 



df 



(1 + ecos /)2 
(21) 



piQ) 



3ttQ 



(22) 



(23) 



^GaAf3(l_e2)3' 



Solving for Q, I get 



37r 3 



(24) 



The uncertainty in Q can be conservatively assessed by linearly adding 
the various sources of errors as 



SQ < 6Q\a + 6Q\g + SQ\m + 6Q\p, + 6Q\e, (25) 

^It agrees with the expression of the anomalistic period of a sateUite orbiting an oblate 
planet obtained in (Capderou, 2005): for a direct comparison Q = —MB?J2, where J2 is 
the first even zonal harmonic of the multipolar expansion of the Newtonian part of the 
gravitational potential of the central body of mass M and equatorial radius R. 
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with 



SQ\a< 



SQ\g< 



Sa, 



I\ / aM3(l-e2)» 
67r V G 



5G, 



SQ\m< |^VG'Ma(l-e2)3-|a2(l-e2)3/2 



(5M, 



^GaM3(l-e2)3 
Stt 



(5Q|e< 



(26) 
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